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In this article, we review the research on the dynamics of quantized vortices in superfluid helium 
and rotating Bose-Einstein condensates with emphasis on the recent research done by our group. 
A quantized vortex is a topological defect that arises from the order parameter in Bose-Einstein 
condensation in which frictionless superfluid flows with quantized circulation around each vortex. 

A quantized vortex was both predicted and discovered first in superfluid *He which was the 
first example of a Bose-Einstein condensate. Quantized vortices have been thoroughly studied in 
superfiuid *He; one of the principal problems was the superfiuid turbulent state consisting of a 
tangle of quantized vortices in thermal counterfiow. More recently, the interest has shifted to the 
nature of superfiuid turbulence, apart from the case of counterfiow. After briefly reviewing the 
earlier research and describing the current problems, we focus our review on superfluid turbulence 
and vortex filament dynamics. One of the important problems is how superfluid turbulence relates 
to classical turbulence. Superfluid turbulence was recently shown to have an energy spectrum 
consistent with the Kolmogorov law, which is an important statistical law in fully developed classical 
turbulence. We also describe the diffusion of an inhomogeneous vortex tangle with relation to the 
observed decay of vortices at very low temperatures where the normal fluid component is so negligible 
that the usual mutual friction does not work as a decay mechanism. In connection to the above 
discussion of groups of vortices, we describe the vortex states that appear in a rotating channel 
with counterflow. Rotational effects cause the vortices to form ordered arrays, whereas counterfiow 
effects tend to cause disordered vortex tangles; the competition of these two effects makes a new 
state of "a polarized vortex tangle" which opens up superfluid phase diagrams as a new area of 
study. 

In the speciflc field of atomic-gas Bose-Einstein condensation, we discuss recent numerical analysis 
of the Gross-Pitaevskii equation that describes the structure and dynamics of the order parame- 
ter. Consistent with observations, the simulated condensate starts an elliptic oscillation after the 
rotation is turned on, which induces the surface-mode excitations. The vortices develop from these 
surface excitations and then enter the bulk condensate, eventually forming a vortex lattice. When 
a condensate is held in a quadratic-plus-quartic combined potential, the fast rotation makes "a 
giant vortex" in which most vortices are absorbed into a central hole around which the superfiow 
circulates. 

PACS numbers: 



I. INTRODUCTION 

In Bose-Einstein condensations (BECs), a macroscopic 
fraction of bosons constitute a condensate wave function 
from which quantized vortices can form. A quantized 
vortex is a topological defect characteristic of a Bose- 
Einstein condensate that has a quantized circulation of 
superfiow and a thin vortex core of the order of the mi- 
croscopic coherence length. 

Quantized vortices were first predicted theoretically by 
Onsager and Feynman 0. Since their discovery in 
superflu id "^He 0, research on quantized vortices has 
become one of the main subjects in low temperature 
physics . On the other hand, the achievement of Bose- 
Einstein condensation in dilute atomic gases in 1995 
has opened up a new research field, in which quantized 
vortices can be well controlled and visualized using new 
optical techniques. The study of quantized vortices has 
an important role in understanding the physics in both 
the traditional field of superfluid helium and the novel 
field of atomic-gas Bose-Einstein condensates (BECs). 



Being motivated by the recent developments in these 
fields, our group has studied theoretically and numer- 
ically the dynamics of quantized vortices in superfluid 
■*He and atomic BECs. The purpose of this article is to 
review the physics of quantized vortices in both systems 
with emphasis on the recent activities in our group j^. 

Section ITU summarizes the properties of quantized vor- 
tices and their differences from vortices in a classical fluid. 
Section IIIII is concerned with the recent research of su- 
perfluid turbulence. Most of the early studies of super- 
fluid turbulence were limited to the case of thermal coun- 
terflow, which has no analogy with classical turbulence. 
However, we also discuss recent experimental studies that 
help clarify the similarity between superfluid turbulence 
and classical turbulence. After describing the equation 
of motion of quantized vortex fllaments, we discuss the 
energy spectrum of a vortex tangle in superfluid turbu- 
lence, the diffusion of an inhomogeneous vortex tangle, 
and a vortex tangle polarized by rotating the vessel. The 
vortex dynamics and vortex lattice formation of rotat- 
ing BECs are described in Sec II VI Here we discuss the 
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theoretical and numerical analysis of the problem. Re- 
cent experimental and numerical results show that el- 
liptic deformation of the condensate and surface-mode 
excitations cause a vortex lattice to form. Section is 
devoted to conclusions. 

II. WHAT IS A QUANTIZED VORTEX? 

Below a critical temperature in an ideal Bose g as, a 
finite fraction of the particles occupies the same single- 
particle ground state and forms a BEG. When the par- 
ticles have mutual interaction, single-particle states are 
no longer meaningful. However a condensate wave func- 
tion ^'(r, t) is still defined as the ensemble average of the 
quantum amplitude for removing a particle at position r 
from the condensate. When a BEG is held in an external 
potential V{r), the dynamics of ^'(r,t) is described by 
the Gross-Pitaevskii(GP) equation 

= (-£ V' + + 9mr, t), 

(1) 

where g = ATrh^m/a represents the strength of interac- 
tion characterized by the s-wave scattering length a and 
m is the mass of each particle. Writing "if = \'^\ exp{i6), 
the squared amplitude | ^' ^ is the condensate density and 
the gradient of the phase 9 gives the superfluid velocity 
Vs = {h/m)W9, which is a frictionless flow of the conden- 
sate. As a result, the circulation of Vg around a closed 
path C in the fluid is quantized as 

/ ds-Vs = — (f ds-Ve = nn {n = 0, ±1, ±2, . . .), (2) 
Jc Jc 

with the quantum of circulation k ~ h/m. Such a vortex 
with a quantized circulation is called a quantized vortex; 
it occurs in superfluid '^He, superfluid "^He, and atomic 
BEGs. 

A quantized vortex is different from a vortex in a clas- 
sical viscous fluid. First, the circulation is quantized, 
which is contrary to the classical vortex that can have 
any value of circulation. Second, a quantized vortex is 
a vortex of inviscid superflow. Thus, it cannot decay by 
the viscous diffusion of vorticity that occurs in the classi- 
cal system. A quantized vortex can decay by shortening 
the length of the core through mutual friction with the 
normal fluid, by breaking into smaller and smaller vor- 
tex rings through reconnections and finally changing to 
some elementary excitations, and by transferring energy 
to smaller length scales through a Kelvin wave cascade, 
followed by acoustic emission. Third, the core of a quan- 
tized vortex is very thin, being the order of the coherence 
length defined by ^ = h/ {^2mg\'^\), which is submicron 
in an atomic-gas BEG and only a few angstroms in su- 
perfluid ''He. Because the vortex core is very thin and 
does not decay by diffusion, it is easy identify the posi- 
tion of a quantized vortex in the fluid. These properties 
make a quantized vortex more stable and definite than a 
classical vortex. 



III. DYNAMICS OF QUANTIZED VORTICES 
IN SUPERFLUID HELIUM 



A. EARLY RESEARCH IN SUPERFLUID 
TURBULENCE AND VORTEX DYNAMICS 

Below the A temperature at about 2.2 K at vapor pres- 
sure, liquid ''He enters the superfluid state called helium 
II. Although the condensate is depressed by interatomic 
interactions, this transition is microscopically caused by 
Bose-Einstein condensation. Helium II behaves like an 
irrotational ideal fluid, and its characteristic phenomena 
can be explained well by the two- fluid model The 
two fluid model states that superfluid ''He is a mixture 
of inviscid superfluid and viscous normal fluid with the 
mixing ratio depending on temperature; the superfluid 
corresponds to the BEG 8] and the normal fluid has all 
of the thermal excitations of the system. 

Most early experimental studies focused on thermal 
counterflow. When helium II is confined in a channel 
closed at one end with a heater, superfluid enters from 
the open end and flows toward the heater. On reach- 
ing the heater, normal fluid is created, which then flows 
back toward the open end. This situation is called ther- 
mal counterflow. It is laminar at low relative velocities, 
but when the relative velocity exceeds some critical veloc- 
ity, the superflow becomes turbulent and has dissipation. 
The concept of superfluid turbulence was introduced by 
Feynman 0] who proposed that the superfluid turbulent 
state consists of a disordered set of quantized vortices, 
which is called a vortex tangle, that dissipates via mu- 
tual friction between the vortex cores and the normal 
flow. The mutual friction was introduced by Gorter and 
Mellink during their study of heat conduction through 
helium II in which they measured the relationship be- 
tween temperature gradient and heat-current density |^ . 
Later, independently of Feynman's proposal. Hall and Vi- 
nen started the pioneering experimental works by making 
the simultaneous measurement of the temperature gradi- 
ent and the second sound attenuation across counterflow. 
They discovered extra second sound attenuation and un- 
derstood it in terms of the two fluid model with a gen- 
eralized mutual friction [l^. After Feynman's proposal, 
Hall and Vinen confirmed that mutual friction occurs by 
observing the anisotropy of second sound attenuation in 
uniformly rotating helium II Later, Vinen investi- 

gated in detail the temperature gradient and the second 
sound attenuation across thermal counterflow to under- 
stand how the mutual friction changed when the heat 
current was suddenly changed. He found that the mu- 
tual friction came from rotational motion of superfluid, 
eventually establishing the picture of superfluid turbulent 
state consistent with Feynman's picture |l2l |. The quan- 
tum of circulation was also measured by Vinen when he 
investigated the oscillation of a vibrating wire in rotating 
helium II There were many experimental studies of 
superfluid turbulence after Vinen's papers ;i33, chiefly on 
thermal counterflow. 
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The nonlinear and nonlocal dynamics of vortices long 
delayed progress in further microscopic understanding of 
the vortex tangle. But Schwarz overcame these difhcul- 
ties 01 . His most important contribution was his di- 
rect numerical simulation of vortex dynamics connected 
with dynamical scaling analysis, thus enabled others to 
calculate such physical quantities as the vortex line den- 
sity, various anisotropic parameters, and the mutual fric- 
tion force 01 . The observable quantities obtained with 
Schwarz's theory agreed well with the experimental re- 
sults of the steady state of the vortex tangle. This re- 
search field pioneered by Schwarz has generated many 
new areas of study in vortex dy namics, such as the 
flow properties in channel s 11^ iWL iVA [T^ , sideband in- 
stability of Kelvin waves Il9j, vortex array in rotating 
SUP er fluid |2 0(1 ■ and vortex pinning [2ll.l2^. 



B. RECENT STUDIES ON SUPERFLUID 
TURBULENCE 

Although our understanding of superfluid turbulence 
and vortex dynamics has made significant gains, the re- 
lation between superfluid turbulence and classical turbu- 
lence remains a major unsolved problem [23L l2^ . Prob- 
ably, this situation has arisen because most studies of 
superfluid turbulence have been devoted to thermal coun- 
terflow which has no analogy with classical turbulence. 

We now discuss several relevant features of classical 
turbulence. In particular, we describe the finding that 
superfluid turbulence can give a clearer example of the 
inertial range than does classical turbulence. It has been 
known for over a hundred years that a classical viscous 
fluid becomes turbulent when the Reynolds number ex- 
ceeds about 1000 25]. The Reynolds number is given by 
the ratio of the nonlinear inertial term to the viscous term 
in the Navier-Stokes equation; therefore, a large Reynolds 
number means that viscosity is less important compared 
with the inertial term. Apart from some boundary ef- 
fects, at high Reynolds numbers, fully-developed, ho- 
mogeneous turbulence follows some universal statistical 
laws. It turns out to be more helpful to consider en- 
ergy and velocity in wave number k space based on the 
spatial Fourier transform v{k) of the velocity field v{r). 
The important statistical law is represented in the energy 
spectrum E{k), where fc = |A;| characterizes a length scale 
and E(k)dk is the average turbulent energy per unit 
mass in the range of the wave numbers dk. As an exam- 
ple, consider stationary grid turbulence in which a fiow 
passing through a grid with the mesh size D is close to 
being homogeneously turbulent far behind the grid. The 
behavior of E{k) is classified into three ranges depending 
on the wave number. In the first energy-containing range, 
the energy is injected into the velocity field at the wave 
number k ^ D~^, by a rate e per unit mass. This en- 
ergy flows into the second inertial range with higher wave 
number; this inertial range is characteristic of fully devel- 
oped turbulence. Here the energy, not being dissipated. 



just flows toward the higher k region by the energy flux e. 
The inertial range has no characteristic scale because it 
is sustained by the Richardson cascade process in which 
large eddies are continuously broken up self-similarly to 
smaller ones by a nonlinear interaction. Since the inertial 
range should be characterized only by the energy flux e 
and the wave number fc, the energy spectrum follows the 
well known Kolmogorov law 

E{k) ^ CKe^^^k-'/^ (3) 

Here, the Kolmogorov constant Ck is found to be about 
1.5 experimentally. The energy flowing from the inertial 
range into the third energy-dissipative range is flnally dis- 
sipated by the dissipation rate e at about the Kolmogorov 
wave number k^: = (^jv^)^!^ with the kinematic viscos- 
ity V. The important point is that the inertial range 
has a universal spectrum that is independent of the spe- 
ciflc dissipative mechanism that operates in the energy- 
dissipative range. It is for this reason that superfluid 
turbulence is expected to follow the Kolmogorov law. 

Recent experimental research studied superfluid tur- 
bulence not in thermal counterflow, thus flnding support 
for the Kolmogorov law. Maurer and Tabeling measured 
local pressure fluctuations in helium flows driven by two 
counter-rotating disks in a range of temperature between 
1.4 and 2.3 K and obtained the Kolmogorov spectrum 
above and below T\ |2^. A group at the University of 
Oregon reported in a series of papers ^7, 28, 22., the 
observed attenuation of second sound behind a grid that 
moved steadily through helium II at temperatures above 
1 K. Among these works, Stalp et al. measured the decay 
of grid turbulence in helium II and showed that the exper- 
imental results were consistent with a classical model of 
energy spectrum that included the Kolmogorov law . 
Then, Vinen analyzed the similarity between superfluid 
turbulence and classical turbulence and showed the 
importance of length scales for understanding the energy 
of the velocity field. For example, although superfiuid 
turbulence is made of a tangle of quantized vortices, the 
situation depends on whether the length scale is larger 
or smaller than the vortex line spacing £. Furthermore, 
at relatively high temperatures, the normal fluid and the 
superfluid are coupled by mutual friction at all relevant 
scales larger than £, and hence the system follows the 
Kolmogorov law at these scales. This theory was shown 
to be consistent with the previous experimental results 
by Stalp et al., thus forming a consistent picture of super- 
fluid turbulence with an appreciable component of nor- 
mal fluid. 

Then appears a very important question; how is the en- 
ergy spectrum of superfluid turbulence at very low tem- 
peratures where the normal fluid component is negligi- 
ble? This is significant in the following reasons. First, su- 
perflow guarantees that any dissipation could work only 
at some very small scale because of its inviscosity. This 
means the presence of the definite inertial range in which 
the energy is transferred from large scale to smaller scales 
by the Richardson cascade process. The resulting small 
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vortices whose size becomes close to the coherence length 
would be unable to keep its vortex nature to change 
into some elementary excitations, or dissipate by acoustic 
emission ^2^, but the inertial range should be indepen- 
dent of such dissipative mechanisms. Second, superfluid 
turbulence is made of a tangle of quantized vortices. The 
inertial range of a classical fluid is believed to be sus- 
tained by the Richardson cascade, while the identifica- 
tion of each vortex in a turbulence is rather obscure. On 
the other hand, quantized vortices are definite and stable 
as described in Sec. II, so that the physical picture of the 
inertial range and the Richardson cascade could be much 
clearer than that in a classical turbulence. Therefore su- 
perfluid turbulence without the normal fluid component 
may give a typical and simple prototype of turbulence. 

There are no experimental studies of the energy spec- 
trum of superfluid turbulence at very low temperatures. 
Nevertheless, the two numerical works on this topic in- 
dicate that the Kolmogorov law also applies to super- 
fluid turbulence at very low temperatures. Starting from 
a flow that mimicked a Taylor-Green vortex, Nore et 
al. studied the decaying turbulence by using the GP 
equation|3l|. The energy spectrum succeeded in showing 
a transient Kolmogorov form over a range of wave num- 
bers less than £^^, but the acoustic emission is closely 
connected with the vortex dynamics and the situation is 
complicated. Our group studied the energy spectrum un- 
der the vortex filament formulation and found that the 
Kolmogorov law should apply to superfluids at very low 
temperatures 32]. This formulation and the results will 
be described in the following subsections. 



C. DYNAMICS OF QUANTIZED VORTEX 
FILAMENTS 

All numerical calculations in this chapter were done 
under the vortex fllamcnt formulation. This subsection 
describes the equations of motion of quantized vortex fll- 
aments and the method of the numerical simulation. A 
vortex fllament is an idealized model of rotational flow 
in which the vorticity is confined to a small core region 
around a one-dimensional line embedded in the three- 
dimensional flow ■ In the field of classical fluid dy- 
namics, vortex fllaments have been a useful model since 
the early 1980s for understanding the geometry and dy- 
namics of a flow. However, in very recent years interest 
in this model has decreased because it is not capable of 
describing either the complicated turbulent structures of 
real systems or solutions of the Navier-Stokes equation 
[25I I . But contrary to the situation in a classical fluid, vor- 
tex filaments are realistic in superfluids such as helium 
II, because here the cores are very thin and relatively 
stable. A quantized vortex filament may be likened to 
a "skeleton" of a vortex because the quantum mechan- 
ical constraint means that each vortex has the ground 
state circulation and a simple potential fiow around its 
core. Our present main interest is to understand if a tan- 



gle of such quantized vortex filaments can still create the 
statistical law such as the Kolmogorov spectrum. This 
subsection is devoted to the description of the vortex fil- 
ament formulation Hi H, followed by a discussion of 
the energy spectrum in the next subsection. 

The vortex filament formulation represents a quantized 
vortex as a filament passing through the fiuid and having 
a definite direction corresponding to its vorticity. Except 
for the thin core region, the superflow velocity field has 
a classically well-defined meaning and can be described 
by ideal fluid dynamics. The velocity at a point r due to 
a fllament is given by the Biot-Savart expression 
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(4) 



where k is the quantum of circulation. The filament is 
represented by the parametric form s — with the 

one-dimensional coordinate ^ along the filament. The 
vector Si refers to a point on the filament and the inte- 
gration is taken along the filament. Helmholtz's theorem 
for a perfect fluid states that the vortex moves with the 
superfluid velocity. Attempting to calculate the velocity 
Vs at a point r = s on the filament makes the integral 
diverge as Si — > s. To avoid this divergence, we sepa- 
rate the velocity s of the filament at the point s into two 
components 
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\si - r\ 



. . (5) 

The first term is the localized induction field arising from 
a curved line element acting on itself, and ^+ and 
are the lengths of the two adjacent line elements, af- 
ter discritization, that hold the point s between them. 
The prime denotes differentiation with respect to the arc 
length ^. The mutually perpendicular vectors s', s" and 
s' X s" point along the tangent, the principal normal and 
the binormal at the point s, respectively, and their mag- 
nitudes are 1, and with the local radius R of 
curvature. The parameter ao is a cutoff parameter equal 
to the core radius. Thus, the first term represents the 
tendency for the points to move the local point s along 
the binormal direction with a velocity inversely propor- 
tional to R. The second term represents the nonlocal field 
obtained by integrating the integral of Eq. along the 
rest of the filament. 

The approximation that neglects the nonlocal terms 
and replacs Eq. (O by s = (3s' x s" is called the localized 
induction approximation (LIA). Here the coefficient (3 is 
defined hy [3 = (k/Att) In (c(i?) /ao), where c is a constant 
of order 1 and is replaced by the mean radius 

of curvature (R) along the length of the filament. Most 
of Schwarz's numerical studies on vortex tangles used the 
LIA because this approximation can greatly reduce com- 
putation times. Although the method is effective for the 
analysis of dense tangles (due to cancellations between 
nonlocal contributions), it does not include the intervor- 
tex interaction properly. Because our present problems 
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described in the following subsections need to take much 
account of the interaction, our numerical simulations are 
done by not the LIA but the fully Biot-Savart law of Eq. 

m ■ 

A better understanding of vortices in a real system re- 
sults when one includes the boundaries in the analyses. 
For this, the boundary-induced velocity field Vg^i, is added 
to Vg so that the superflow can satisfy the boundary con- 
dition of an inviscid flow. If the boundaries are specular 
plane surfaces, v^^b is just the field due to an image vortex 
made by reflecting the vortex into the plane and reversing 
its direction of vorticity. To allow for another, presently 
unspecified, applied field, we include Vs,a- Hence, the to- 
tal velocity so of the vortex filament without dissipation 
is 



So 



— s'x s"ln 

4tt 



.1/4 
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At finite temperatures the mutual friction between the 
vortex core and the normal flow Vn is important. Includ- 
ing this term, the velocity of s is given by 

s = So + as' X {vn - So) - a's' x [s' x (v„ - Sq)], (7) 

where a and a' are the temperature-dependent friction 
coefficients, and Sq is calculated from Eq.ljHJ. 

The method used in the numerical simulation is simi- 
lar to that of Schwarz jlj| and described in detail in our 
paper [33 |. A vortex filament is represented by a single 
string of points with a distance A^. The vortex con- 
figuration at a given time determines the velocity field 
in the fluid, thus moving the vortex filaments according 
to Eqs. l|njl and (TJ. Both local and nonlocal terms are 
represented by means of line elements connecting two ad- 
jacent points. As the vortex configuration develops and, 
particularly, when two vortices approach each other, the 
distance between neighboring points can change. Then 
it is necessary to add or remove points properly to retain 
sufiicient local resolution. Through the cascade process, 
a large vortex can break up through many reconnections, 
eventually becoming a vortex that is smaller than the 
space resolution A^. The simulations cannot follow the 
dynamics beyond this point, so these vortices are elimi- 
nated numerically. 

It is important to properly include vortex reconnection 
when simulating vortex dynamics. A numerical study of 
a classical fluid showed that the close interaction of two 
vortices leads to their reconnection, chiefly because of 
the viscous diffusion of the vorticity [ssj . Schwarz as- 
sumed that two vortex fllaments reconnect when they 
get close within a critical distance, and showed that the 
statistical quantities such as vortex line density were not 
sensitive to how to make reconnections. Even after the 
Schwarz's works, it was still unclear whether quantized 
vortices can actually reconnect or not. However Koplik 
and Levine solved directly the GP equation to show the 



two close quantized vortices reconnected even in an in- 
viscid superfluid [s^. More recent simulations showed 
that reconnections were accompanied by acoustic emis- 
sions [s^, nil . As the vortex filament formulation cannot 
represent the reconnection process itself, we assume that 
two vortices reconnect when they approach within the 
space resolution Af . The detail of this procedure is dis- 
cussed in our paper [s^. 



D. ENERGY SPECTRUM OF SUPERFLUID 
TURBULENCE 

To understand the energy spectrum and determine 
whether or not superfluids follow the Kolmogorov law 
at very low temperatures, one must calculate the energy 
spectrum. The energy spectrum is originally calculated 
by the Fourier transform of the fluid velocity v{r). The 
superfluid velocity Vs{r) is determined by the configura- 
tion of quantized vortices in our vortex filament formu- 
lation. Therefore, one can calculate the energy spectrum 
directly from the configuration of vortices. This is a big 
advantage of this formulation because it makes the nu- 
merical calculation of the spectrum less time-consuming 
than the method via v{r). Using the Fourier transform 
Vs{k) = {2tt)~^ J drVs{r) ex-p{—ik-r) and Parseval's the- 
orem J dk\vs{k)\'^ = (27r)~^ J dr\vs{r)\'^ , the kinetic en- 
ergy of the superfluid velocity per unit mass is 



E 



1 



dr\vs{r)\'' 



dk\v,{k)[^ 



(8) 



The vorticity u){r) = lotv s{r) is represented in Fourier 
space as Vs{k) = ik x u){k)/\k\'^, so that we have 
E = ((27r)V2)/dfc|u;(fc)|V|fe|2. The vorticity a;(r) = 
K J (i^s'(^)^(s(^) — r) in the vortex filament formulation 
is rewritten as uj{k) — (K/(27r)^) / s'(^) exp(— is(^) • 
fe). Using the definition of the energy spectrum E{k) 
from E — dkE{k), these relations yield 

xexpHfc.(s(ei)-s(e2))), (9) 

where dflk — sm9kd6kd4>k is the volume element in 
spherical coordinates. This formula connects the energy 
spectrum directly with the vortex conflguration. 

Starting from the Taylor-Green vortex and following 
the vortex motion without the mutual friction, we ob- 
tained a roughly homogeneous and isotropic vortex tan- 
gle (Fig. 1) 32]. This is a decaying turbulence, being dis- 
sipated by the cutoff of the smallest vortices whose size 
is comparable to the space resolution A^ = 1.83 x 10~^ 
cm. At flrst, the energy spectrum has a large peak at 
the largest scale where the energy is concentrated, but 
the spectrum changes as the vortices become homoge- 
neous and isotropic. The time dependence of the energy 
dissipation rate e shows that de/dt becomes small after 
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(c) (d) 

FIG. 1: Time development of the vortex tangle ai t — sec 
(a), 30.0 sec (b), 50.0 sec (c) and 70.0 sec (d). The vortices 
are confined in a cube of size D = 1 cm and the calculation 
is made under the fully Biot-Savart law (Araki, Tsubota and 
Nemirovskii, Phys. Rev. Lett. 89,145301-2, 2002, reproduced 
w ith permission. Copyright(2002) by the American Physical 
Society) . 



about 70 sec and artifacts of the initial state disappear 
[32]. Similarly, the isotropic parameters introduced by 
Schwarz J^J indicate a nearly isotropic vortex tangle be- 
comes after 70 sec. Figure 2 shows that the energy spec- 
trum of the vortex tangle at 70 sec agrees quantitatively 
with the Kolmogorov spectrum in the small k region. The 
dissipative mechanism due to the cutoff works only at the 
largest wave number k ^ 27r/A^ = 343 cm~^. However 
the energy spectrum at small k region is determined by 
that dissipation rate. By monitoring the development 
of the vortex size distribution, such decay of a tangle is 
found to be sustained by the Richardson cascade process 
[33| . These results support the classical picture of the 
inertial range in superfluid turbulence at very low tem- 
peratures. 

The power of the spectrum changes from -5/3 to -1 
at about k ^ 2Tr/£. The k^^ spectrum comes from the 
contribution of the velocity near each single vortex. Our 
simulation of this k~^ regime "id] indicates that the en- 
ergy cascade process is due to the Kelvin wave cascade 
process. The simulations showed that when there is con- 
tinuous excitation of small k Kelvin waves along a single 
vortex and if a sink removes the energy at large wave 
numbers, then the nonlinear coupling between different 
modes leads to a net flow of energy from small to large 
wave numbers. This results in a simple steady spectrum 
of Kelvin waves that is insensitive to the strength and 
frequency of the excited drive. 




k [1/cm] 

FIG. 2: Comparison of the energy spectrum (solid line) at 
t = 70 sec with the Kolmogorov law E{k) = Ce^'^^fc"^/^ 
(dotted line) with C = = 1 and e = 1.287 x lO"*^ 

cm^/sec^ law (Araki, Tsubota and Nemirovskii, Phys. Rev. 
Lett. 89, 145301-3, 2002, reproduced with permission. Copy- 
right (2002) by the American Physical Society). 

A specific picture of vortex tangles at very low tem- 
peratures has emerged from these studies. In the low k 
regime, the vortex tangle has the inertial range in which 
the Kolmogorov spectrum is realized and the Richard- 
son cascade process describes the energy flow to larger 
wave numbers; in this regime, superfluid turbulence mim- 
ics classical turbulence. In the regime with k > 2Tr/£, 
the Kelvin wave cascade process becomes relevant to the 
transfer of energy to the high wave numbers above which 
some other mechanism dissipates the energy; this cascade 
process may be obscure in classical turbulence because 
the wave numbers become comparable to the vortex core 
sizes, but it is clearly seen in superfluid turbulence at 
very low temperatures. 

E. DIFFUSION OF AN INHOMOGENEOUS 
VORTEX TANGLE 

At flnite temperatures, vortex tangles decay due to 
mutual friction. We ask here if vortex tangles decay at 
very low temperatures, where the normal fluid compo- 
nent is so negligible that the mutual friction does not 
work, and ask what mechanism removes energy from 
the tangle if it does decay. Experimentally, these ques- 
tions are not easy to answer because counterflow and 
second sound techniques cannot work at such low tem- 
peratures. Nevertheless, Davis et al. created vortices in 
superfluid ''He using a vibrating grid, and observed their 
temperature-independent decay below T ^ 70 mK by 
the trapping of negative ions at vortices ^3 ■ The energy 
cascade process described in the previous subsection may 
be closely related to this decay, although an understand- 
ing of the behavior of an inhomogeneous vortex tangle is 
likely needed. Thus, we studied numerically the diffusion 
of a localized tangle 14l|. We flrst prepared a localized 
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vortex tangle. Then we followed its decay and diffusion 
under the fully Biot-Savart law in the absence of mutual 
friction. 

The decay of a homogeneous vortex tangle is described 
by Vinen's equation dL/dt — —X2{i^/2tt)L'^ ,4^]. Here 
L = (1/A) j d^ is the vortex line density, where ^ is the 
one-dimensional coordinate along the vortex filaments 
and the integral is taken along all vortices in the sample 
volume A. The parameter X2 depends on temperature, 
being about 0.3 at zero temperature 0|. A simple gen- 
eralization of this equation for an inhomogeneous system 
may be obtained by adding a diffusion term, so that we 
have the inhomogeneous Vinen's equaion 

^^-X2^L{r,tf + DW'L{r,t), (10) 

dt ZTT 

where Z? is a diffusion constant. A comparison between 
the numerical solution of Eq. H10|l and the above nu- 
merical simulation shows that D — (0.1 ± 0.05)k. Hence 
the diffusion of a tangle is relatively small; similar results 
were obtained by Barenghi and Samuels 43|. Even if a 
tangle is created by a vibrating grid, the diffusion is small 
and the tangle can remain localized in the neighborhood 
of the grid as far as there is no motion on a length scale 
larger than the line spacing. If motion on a larger scale 
is generated, then the diffusion can be enhanced. 

F. ROTATING SUPERFLUID TURBULENCE 

Rotating superfluid turbulence is interesting because 
the rotation can cause vortex ordering and thus lead to 
various vortex phases. This problem is different from 
that of superfluid turbulence at zero temperature, but 
closely related with the motivation becoming very im- 
portant recently. Most configurations of quantized vor- 
tices which have been investigated in superfluid helium 
can be grouped into two types: ordered vortex arrays 
and disordered vortex tangles, and the former part of 
this article has been devoted to the disordered vortex 
tangle. For example, ordered arrays of vortices occur 
when superfluid helium is rotated with angular veloc- 
ity fl exceeding a certain critical value Q. The result- 
ing quantized vortices are aligned along the rotation axis 
and form an array whose areal number density is 2f2/K, 
in agreement with a rule discovered by Feynman. Thus, 
an important question arose: what happens if vortices 
are created by both rotation and counterflow along the 
rotational axis? Only one experiment seems to address 
this issue 44]. In the experiment, Swanson, Barenghi, 
and Donnelly mounted a counterflow channel on a rotat- 
ing cryostat, thus being able to create vortices by inde- 
pendent combination of rotation and counterflow. The 
vortex line density L was determined from the second- 
sound attenuation across the channel. They found two 
critical velocities with Vd < Vc2- At the lowest rotation 
velocities, when the counterflow velocity Vns is less than 
Vci, a vortex array formed with a density in agreement 



with Feynman's rule. This critical value of Vd was much 
less than that with no applied rotation. In addition, the 
measured Vd was consistent with the onset of a vortex 
wave instability discovered experimentally by Cheng et 
al. and explained by Glaberson et al. as be- 

ing caused by helical Kelvin waves along the vortex cores 
that are destabilized by the component of counterflow ve- 
locity along the vortices. This instability is usually called 
the Donnelly-Glaberson (DG) instability. The nature of 
the flow within this unstable regime and the cause of the 
second critical velocity Vc2 is not clear. 

To better understand these flow regimes, we simulated 
the original experiments of Swanson et al. !43| , in which 
the superfluid flowed in a rotating cube, and found that 
the ordering of the vortex array was sensitive to the rota- 
tion velocity. In a rotating cube, the equation of motion 
of vortices is modifled by two effects. The first effect 
is the force acting upon the vortex due to the rotation. 
According to Helmholtz's theorem, the generalized force 
acting upon the vortex is balanced by the Magnus force 
as Psk{s' X So) = SF' /Ss, where ps is the superfluid den- 
sity, F' — F — ft ■ M is the free energy of a system in the 
frame rotating with the angular velocity Jl, and M is 
the angular momentum. Taking the vector product with 
s', we obtained the velocity Sq. The first term F due to 
the kinetic energy of vortices gives the Biot -Savart law of 
Eq. ©. The second term • M leads to the velocity Sj-ot 
of the vortex caused by the rotation. The second effect 
is the superflow induced by the rotating square cross- 
section. For a perfect fluid, we can use the analytical so- 
lution of the velocity Vcuh inside a cube of size D rotating 
with angular velocity Jl = Viz ^^l- By including both 
effects, we obtained the velocity Sq in a rotating cube as 
the sum of v^uh and Sj-ot, the right-hand side of Eq. jSJ. 
In the numerical simulation, the sample was a cube of 
size Z? = 1.0 cm. We used periodic boundary conditions 
along the rotating axis and rigid boundary conditions at 
the side- walls. Also, the counterflow was applied along 
the z-axis, the normal fluid was assumed to be at rest 
in the rotating frame, and, to make comparison with the 
experiment |[4Jj , we did the calculation for a temperature 
T — 1.6 K. Hence, the calculation included mutual fric- 
tion. Two quantities were introduced to characterize the 
vortex tangle. One is the vortex line density L, and the 
other is the polarization of the tangle which is defined as 
(s^) = (1/AL) / d^s' ■ z. Thus (s^) is unity for a perfect 
vortex array and zero for a randomly oriented tangle. 

The numerical simulation began from a rotating vortex 
array with small random perturbations added to the vor- 
tices, then a counterflow velocity was applied. Figure 3 
shows how 33 vortices, initially parallel to each other in a 
stable configuration at = 4.98 x 10~^ rad/sec, develop 
into a vortex tangle under T4s = 0.08 cm/sec. The par- 
allel vortices become unstable and Kelvin waves grow by 
the DG instability. When the amplitude of the Kelvin 
waves becomes comparable to the average vortex sepa- 
ration, vortex reconnections take place and more vortex 
loops are created. These loops disturb the initial vortex 
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FIG. 4: Dependence of L and (s^) on v^s in the saturated 



state for (1 = (white square), f2 = 9.97 x 10 rad/sec (black 
square), SI = 2.99 x 10"^ rad/sec (circle) and Q. = 4.98 x 10"^ 
rad/sec (triangle). The right figure lacks the data for £7 = 
(Tsubota, Araki and Barenghi, Phys. Rev. Lett. 90, 145301- 
3, 2003, reproduced with permission. Copyright(2003) by the 
American Physical Society). 



IV. DYNAMICS IN ROTATING 
BOSE-EINSTEIN CONDENSATE 



FIG. 3: Development from a vortex array to a polarized tangle 
at T = 1.6 K,n = 4.98 X 10"^ rad/sec and Vns = 0.08 cm/sec. 
Computed vortex tangle at the following times: (a): t=0 sec; 
(b): 12sec; (c): 28 sec; (d): 160 sec (Tsubota, Araki and 
Barenghi, Phys. Rev. Lett. 90, 145301-3, 2003, reproduced 
with permission. Copyright(2003) by the American Physical 
Society) . 



array, leading to an apparently random vortex tangle. 

The results show that rotation tends to order the vor- 
tex tangle over that which would occur without rotation 
[T^ 13 . The values of L and (s^) show that the vor- 
tex tangle evolves into a statistically saturated steady 
state. A vortex tangle created with only counterflow is 
sustained by the competition between the driving coun- 
terflow and the mutual friction 14], whereas in our case 
the steady state is achieved by balancing the effects of ro- 
tation, counterflow, and mutual friction. Figure 4 shows 
the dependence of L and (s^) on V^^^ and in each satu- 
rated state. It is evident that the polarization decreases 
with increasing counterflow velocity and increases with 
increasing rotation. This shows the competition between 
ordering caused by rotation and disordering caused by 
counterflow. 

Recently Finne et al. used NMR to study turbulence 
of the B phase of superfluid ■^He and observed a sharp 
velocity-independent transition at a critical temperature 
between two regimes Regular behavior occurred at 
high temperature and turbulence occurred at low tem- 
peratures. They also found that the experimental results 
were consistent with numerical vortex dynamics simu- 
lations which was done by the method described here. 
This turbulence comes from intrinsic instability of vor- 
tices depending on temperature, being different from the 
DG instability. 



A. QUANTIZED VORTICES IN 
ALKALI-ATOMIC BOSE-EINSTEIN 
CONDENSATE 

The recent dramatic achievement of Bose-Einstein con- 
densation in trapped alkali-atomic gases at ultra-low 
temperatures has stimulated intense experimental and 
theoretical activity Such atomic-gas Bose-Einstein 
condensates (BECs) differ fundamentally from liquid he- 
lium condensates in several ways. First, the condensates 
of alkali-atomic gases are dilute, having the mean par- 
ticle density n with n\a\'^ ^ 1, so that the interatomic 
interaction can be accurately parametrized in terms of a 
scattering length a. As a result, at low temperatures, the 
GP equation ^ gives an extremely precise description 
of the atomic condensate and their dynamics. This sit- 
uation differs from superfluid ^He, where the relatively 
high density and strong repulsive interactions complicate 
greatly the analytical treatments. Furthermore, because 
of the relatively strong interactions, the condensate frac- 
tion in bulk superfluid *He is only about 10% of the total 
particles, even at zero temperature. In contrast, almost 
all atoms participate in the condensate in an atomic-gas 
BEG. Second, a BEG in bulk helium has uniform den- 
sity, whereas an atomic-gas BEG has nonuniform density 
due to the confinmcnt by a trapping potential. In this 
section, we review the dynamics of quantized vortices in 
an atomic-gas BEG system. 

From the above properties, the physics of quantized 
vortices has novel features. First, since the density is 
dilute, the relatively large coherence length ^ makes it 
possible to directly visualize the quantized vortices by 
using optical techniques. Second, because the order of 
the coherence length ^ is close to the size of the con- 
densate, the vortex dynamics is closely connected with 
the collective motion of the condensate density itself. As 
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described later, vortex nucleation of this system is re- 
lated to an instability in the surface excitations of the 
condensate, which differs from that in superfluid helium 
system where the compressibility of the density is proba- 
bly negligible. Finally, alkali atoms have internal degrees 
of freedom attributed to the hyperfine spin, so that their 
BECs can have multi-component order parameters if the 
spin degrees are available. Multicomponent BECs pro- 
vide us with new possibilities to study unconventional 
vortex states that have been studied extensively in other 
fields of physics such as superfluid '^He, anisotropic su- 
perconductors, and cosmology. 

Superfluids can support dissipationless flow, which is 
closely related to the existence of quantized vortices. 
Hence, the importance of studying a quantized vortex 
was recognized immediately after the first atomic-gas 
BEC was realized, in order to know whether or not the 
system becomes superfluid. The first experimental de- 
tection of a vortex involved a ^^Rb BEC containing two 
different internal hyperfine components . The genera- 
tion of a vortex was achieved by a phase engineering tech- 
nique, in which interconversion between two components 
was controlled spatially and temporally by an external 
couphng field |5lj. They then turned off the coupling, 
thus creating a vortex state consisting of one circulating 
component surrounding a nonrotating core of the other 
component. The ENS group observed the formation of 
single and multiple vortices in a single-component , *^Rb 
elongated cigar-shape condensate |52|- The condensate 
was trapped in a static axisymmetric magnetic trap and 
a nonaxisymmetric attractive dipole potential created by 
a stirring laser beam. This combined potential produces 
a cigar-shaped harmonic trap with a slightly anisotropic 
transverse profile. By rotating the transverse anisotropy 
at a frequency Q, they could observe the formation of a 
vortex above a certain critical value of A lattice con- 
sisting of more vortices appeared when fl was increased 
further. In contrast to the indirect visualization methods 
used in superfluid helium systems, the quantized vortices 
were directly visualized as holes in the transverse density 
profile of the time-of-flight images. Later, groups at MIT 
[53, JILA [s^l, and Oxford IK^ observed vortex lattices; 
the first two groups observed lattices that consisted of 
100 or more vortices. 

There are numerous theoretical studies of quantized 
vortices in atomic-gas BECs (see Ref. H^)- In this 
area, our group focused on the dynamical properties of 
quantized vortices. Based on the numerical simulation of 
the time-dependent CP equation, we revealed the mar- 
velous nonlinear dynamics that a rotating BEC under- 
goes li^llEll^- Ssc. IIVBI we discuss the dynamics 
of vortex lattice formation in a rotating BEC |53, |59|, 
motivated by the experimental observation of Ref. [6Clj| . 
This work clarified for the first time the nonlinear dy- 
namics of vortex lattice formation, in good agreement 
with the experimental results. In Sec. IIV CI we study 
the dynamics of a fast rotating condensate ^58] . By com- 
bining an additional quartic potential with a harmonic 



potential, we find the new vortex states involving "a gi- 
ant vortex" , in which many vortices are absorbed in a 
single density hole. 



B. VORTEX LATTICE FORMATION IN A 
ROTATING BOSE-EINSTEIN CONDENSATE 

Madison et al. observed directly nonlinear dynamical 
phenomena such as vortex nucleation and lattice forma- 
tion in a rotating condensate |63|- By suddenly turning 
on the rotation of the potential, the initially axisymmet- 
ric condensate made a collective quadrupole oscillation 
in which the condensate deformed elliptically. This os- 
cillation continued for a few hundred milliseconds with 
a gradually decreasing amplitude. Then, the axial sym- 
metry of the condensate suddenly reappeared and con- 
currently the vortices entered the condensate from its 
surface, eventually settling into a lattice configuration in 
the bulk. Such a direct observation of vortex nucleation 
and vortex lattice formation has never been done in su- 
perfluid helium. This observation motivated us to study 
theoretically the detailed dynamics of vortex lattice for- 
mation. 

One motivation for studying a rotating, atomic-gas 
BEC is to determine whether the dynamical instability 
|6ll |. the Landau instability J^j, or both are responsible 
for vortex nucleation and lattice formation. The dynami- 
cal instability originates from the imaginary frequency of 
the excitation mode, giving rise to an exponential growth 
of the unstable mode even in the energy-conserving dy- 
namics. The Landau instability occurs when the exci- 
tation spectrum has negative eigenvalues in the rotating 
frame and the system is subject to some energy dissi- 
pative mechanism. These two instabilities usually oc- 
cur in different parameter regimes. Thus, one may ask 
which instability is important for actual vortex nucle- 
ation. In the experiments ^2, 53, 60], the vortices nu- 
cleated most frequently when the rotation frequency O 
was near 0.7oj±, where uj± is the transverse trapping fre- 
quency. The rotating potential in these experiments ex- 
cites mainly the surface mode with the quantum number 
of angular momentum I — 2 (quadrupole mode). When 
the interaction energy of the condensate is much larger 
than the energy of a harmonic potential, the dispersion 
relation for the surface mode is given by uji = Vlujj_ 
[63l |. In the rotating frame with frequency fi, the surface- 
mode frequency is shifted by —l^l. Hence, it is expected 
that the quadrupole mode I — 2 is resonantly excited 
at il/uj± = \/2/2 ~ 0.707. Sinha and Castin showed 
that the rotating condensate without vorticity was dy- 
namically unstable near the quadrupole resonance and 
pro posed that this instability triggered vortex nucleation 
[6ll |. However, our later analysis reveals that both insta- 
bilities are necessary to explain the actual vortex forma- 
tion. 

We resolved these problems by numerically solving the 
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GP equation which is an extension of Eq. [571 159| 

at L Zm J 

Here, the centrifugal term ^QLz — iMl(xdy — ydx) ap- 
pears in a system rotating about the z-axis at a frequency 
n. In this work, we assume translational symmetry along 
the z-axis, thus making the problem two-dimensional. 
Then the normalization of the wave function ip is taken as 
/ drli/jp = n2D with the particle number per unit length 
along the z-axis. The trapping potential has the form 



density 



phase 



V{v) 



—muj' 
2 



[{{l + ex)x^ + {l + ey)y^}, 



(12) 



with the small anisotropy parameters ex and ey {ex ^ ey) 
; this form describes approximately the rotating poten- 
tial used in the ENS experiments [s^ l60| . By intro- 
ducing the scales characterizing the trapping potential 
ttfi = y^h/2mcj± for length and for time, and re- 
placing ip y/n2Bfp/cih, the coupling coefficient in the 
nonlinear term of Eg. Hll|l becomes C — 87ran2D- The 
condition of Ref. l6Qj, corresponding to a = 5.77nm, 
iV = 3 X lO*^, = 11.8 X 2n and A = uj±/uJz = 9.2, gives 
C ~ 500. The term with 7 introduces dissipation, which 
is treated phenomenologically in the GP equation. This 
form of the dissipative equation follows the work of Choi 
et al. j^; they studied the damping of the collective 
oscillation of a BEG and determined the value of 7 to 
be 0.03 by fitting their theoretical results with the MIT 
experiments |6^. Hence, we use 7 = 0.03 here. Because 
this dissipative term is much smaller than other terms 
in the GP equation, a small variation of 7 would not 
change the dynamics qualitatively but can only modify 
the relaxation time scale. This form can be derived from 
the formulation for a finite-temperature BEG developed 
by Zaremba et al. j66j under some approximations [59j |. 
Gardiner et al. obtained the similar but different equa- 
tion by another approach [67j . applying it to the simu- 
lation of vortex lattice formation |68i| . Because the time 
development of Eq. does not conserve the norm of 
the wave function, we had to adjust the chemical poten- 
tial /i at each time step to ensure normalization. 

The simulations started from the stationary solution 
with uniform phase of Eq. Hll|l for a nonrotating trap 
{^x = = 0). The method for turning on a rotating 
drive followed the experimental procedure in Ref. |60j| . 
The rotation with frequency f2 started ai t — 0, and 
the trap anisotoropy e — {(1 4- ex) — (1 + ey)}/{{l + 
^x) + (1 + ey)} was increased rapidly from zero to its fi- 
nal value 0.025 in 20 msec, keeping ey zero. First, we 
discuss the dynamics of the condensate with C — 500 
for n/uj± — 0.7. When the dissipation is free (7 = 0), 
the condensate undergoes a quadrupole oscillation, but 
it shows only the simple periodic oscillation in a sense of 
the Fermi-Pasta- Ulam recurrence [see Fig. 6]. When 
the dissipation is present, this behavior is dramatically 
changed. The left column of Fig. 5 shows the time devel- 
opment of the condensate density \^p{x,y,t)\'^ . Initially, 






FIG. 5: Time development of the density (left column) and 
phase (right column) profile after the trapping potential be- 
gins to rotate suddenly with f2 — 0.7uj± for (a) 67 msec, 
(b) 360 msec, (c) 425 msec and (d) 735 msec. The value of 
the phase varies continuously from (black) to 2tt (white). 
There appear some lines where the phase changes discontin- 
uously from black to white, corresponding to the branch cuts 
between the phase and 2n. The apexes of these lines around 
which the value of the phase rotates continuously from to 
2tt represent the quantized vortices. 



the condensate also undergoes a quadrupole oscillation, 
but the oscillation is damped because of the dissipation. 
After a few hundred milliseconds, the boundary surface 
of the condensate becomes unstable, generating surface 
ripples that propagate along the surface. The excitations 
are likely to occur on the surface whose curvature is low, 
i.e. parallel to the longer axis of the ellipse as shown in 
Fig. 5(b). Then the waves on the surface develop into 
the vortex cores, which enter the condensate. As the vor- 
tices penetrate inside the condensate, the axisymmetry 
of the condensate shape is recovered. As is known in 
the study of rotating superfluid helium [23| , the rotating 
drive pulls vortices toward the rotation axis, whereas re- 
pulsive interactions between vortices tends to push them 
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FIG. 6: Time evolution of the deformation parameter a for 
Q,/lu± — 0.70 with dissipation (bold-solid curve), for n/uj± — 
0.70 without dissipation (dashed-curve) and for fi/ojx — 0.75 
without dissipation (dotted-curve). 



apart; therefore, their competition yields a vortex lattice 
whose vortex density depends on the rotation frequency. 
In this simulation, six vortices form a vortex lattice. Fig- 
ure 6 shows the time development of the deformation 
parameter defined as 



a{t) 



-n- 



(13) 



where (A) means J dxdyip* AiJ; . Madison et al. observed 
that, after a rotation of 17 = 0.7u;_l started suddenly, a 
oscillated during a few hundred milliseconds and then 
reduced abruptly to a value below 0.1 when vortices en- 
tered the condensate from its surface . The evolution 
of a (bold-solid curve in Fig. 6) closely resembles Fig. 
3 of Ref . [gOI , and our result is thus consistent with the 
experimental one. 

The corresponding time development of the phase of 
ip{x,y,t) is shown in the right column of Fig. 5. As 
soon as the rotation starts, the phase field inside the 
condensate takes the form of quadrupolar flow 6{x, y) = 
a(t)xy + const, and just outside the Thomas-Fermi 
boundary appear many proto-vortices; for example. Fig. 
5(a) shows about 20 proto-vortices. When the proto- 
vortices are on the outskirts of the condensate where the 
amplitude |V'(r,t)| is almost negligible, they neither con- 
tribute to the energy nor the angular momentum of the 
system. Because they are invisible in the corresponding 
density profile of F ig. 5, they may be called "ghost vor- 
tices". The ghost vortices move toward the rotation axis, 
but their invasion into the condensate is prevented at the 
boundary surface within which the amplitude grows up. 
However, as the surface ripples are generated, the ghost 
vortices begin to penetrate the condensate. At this point, 
only some defects enter the bulk condensate because their 
further invasion costs energy. For example. Fig. 5(d) 
shows that six vortices enter the condensate and form 
a lattice, while the remaining vortices are repelled and 
escape to the outside of the Thomas-Fermi boundary. 



The penetration of the ghost vortices into the conden- 
sate was accomplished with the help of the surface rip- 
ples, induced by the instabilities in the non- vortex state. 
One possible mechanism for the surface ripple excitation 
is the dynamical instability, which occurs without any 
dissipation. Sinha and Castin did a linear stability anal- 
ysis of a rotating condensate, finding the growth of the 
fluctuation in the resonant range 0.72 < n/uj± < 0.78 
[6l| . To confirm whether this instability leads to the 
vortex lattice formation, we calculated the dynamics for 
il/uj± ~ 0.75, and show the evolution of a in Fig. 6. We 
can see that the perfect recurrence of the oscillation is 
broken. This is due to the dynamical instability, which 
causes the irreversible transfer of the quadrupole mode 
into the higher-energy excitation mode via the nonlinear 
mode coupling. However, although the surface ripples 
were generated, the ghost vortices never penetrated the 
inside of the condensate. Therefore, the dynamical insta- 
bility alone cannot explain the vortex lattice formation. 
Actually, in the experiments, the vortices occurred at 
even lower off- resonant frequencies [s^ Issl Issl l60l | . a find- 
ing that cannot be understood using the dynamical in- 
stability. The other instability is the Landau instability, 
which is applicable to the case with dissipation. The crit- 
ical frequency of this instability can be expressed by the 
Landau criterion flc ~ mm{uji/l) for a rotating BEG [6^ . 
The angular momentum number Ic that yields fie takes 
a value larger than 4 with the parameter used in experi- 
ments 62\. For C = 500, we obtain lc=8 and flc = 0.5uj± 
from uji by numerically solving the Bogoliubov-De Genne 
equation. Thus, our simulation confirmed that this insta- 
bility actually leads to the generation of vortices, where 
the surface ripples always evolved into the vortex cores 
and penetrated into the condensate. 

In the range 0.72 < il/u!± < 0.78 where the dynami- 
cal instability is effective, the surface ripples are excited, 
but dissipation is needed for the vortices to penetrate 
into the condensate. Here the dissipation may originate 
in the thermal component, which should be almost negli- 
gible in the atomic gas experiments at very low temper- 
atures. However, in the dynamical process of a conden- 
sate, there is a possibility that the thermal component 
will be produced under a stron g p erturbation. Hence 
the experimental results of Ref. |52l IssL l60| may be ex- 
plained as follows. Consider a situation in which there is 
almost no dissipation at very low temperatures, namely, 
a case in which the Landau instability cannot occur. In a 
quadrupole resonance, however, the dynamical instability 
causes stochasticity in condensate oscillations, leading to 
the creation of the thermal component. Indeed, Hodby 
et al. has reported that the temperature of the system 
increases from 0.5rc to 0.8Tc when the vortices form [55l| . 
The time spent during this process is determined by the 
growth time of the dynamical instability, which is ex- 
pected to be about 100 msec [gJ- Thus, this oscillation- 
created thermal component makes the dissipation effec- 
tive, so the vortices can penetrate into the condensate 
via the Landau instability. To test this hypothesis we 
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need to use analysis beyond the GP equation and thus 
we leave this issue for future study. In the other range of 
f2, a condensate makes only a stable quadrupole oscilla- 
tion without dissipation. However, our results show that 
vortices may be generated at frequencies above flc when- 
ever finite dissipation occurs. Therefore, if one does an 
experiment in which the temperature is high enough to 
have effective dissipation, one should be able to observe 
the critical frequency given by the Landau criterion. 



C. GIANT VORTEX FORMATION IN A FAST 
ROTATING BOSE-EINSTEIN CONDENSATE 

In this subsection, we consider the dynamics of a ro- 
tating BEC trapped in a harmonic-plus-quartic potential 
[Hsj l. The greater the rotation frequency, the greater the 
number of vortices in the lattice. In this study, we want 
to find out what happens when the size of the vortex 
cores becomes comparable with the vortex separation. In 
such a situation, the vortex lattice is predicted to melt 
via quantum fluctuation, leading to a quantum Hall-like 
state [T^, 0, 01- For a rotating condensate with fre- 
quency 17 in a harmonic trapping potential (l/2)mu;^r^, 
the centrifugal potential cancels the confinement, thus 
preventing a BEC from rotating at Q beyond u!±. Un- 
fortunately, some authors have concluded that the study 
of a dense vortex lattice encounters significant difficulties 
arising from such a softening of the effective trap poten- 
tial [l^JiJ. This difiiculty can be avoided by introducing 
an additional quartic potential, so that the effective po- 
tential for the rotating condensate becomes 



Q=2.5(0, 



V{r) = ^m{ujl 



4 



(14) 



The quartic part of V allows us to increase above 
uj±, then the potential has a "Mexican hat" structure. 
Hence, we can expect the dynamics to be characteristic 
of the quadratic-plus-quartic potential for Q > uj± . As in 
the last subsection, we also consider a two-dimensional 
system subject to rotation = flz by assuming trans- 
lational symmetry along the z axis. In the numerical 
simulation, the rotating drive is suddenly turned on by 
introducing a small anisotropy of the harmonic trap as 
{exx'^+eyy'^)/4:. Finally, the small trap anisotropy (cxjCy) 
is turned off adiabatically to obtain the axisymmctric fi- 
nal steady state. 

Numerical simulations were done for several values of 
at fixed C — 250 and fc = 1. Figure 7 shows the time 
development of the density profiles for CI/luj_ = 2.5 and 
il/uj± = 3.2. As before, the ripples are excited on the 
condensate surface, and some ripples develop into vortex 
cores and penetrate into the condensate. For fl/u;± — 
2.5, however, some penetrating vortices move toward the 
center, merging together to make a density hole at the 
center, around which some other vortices form a circular 
array. The phase profiles of the final steady state are 
shown in Fig. 7(d). Although some phase defects come 



(a) 



(e) 



Q=3.2o), 



(b) 



(f) 



(c) 



olo 




FIG. 7: Time development of the density profile for Q — 
2.50;^ (left column) and = 3. 20)^ (right column) after the 
trapping potential begins to rotate suddenly. The time is (a) 
30 msec, (b) 33 msec, (c) 300 msec for the left column and 
(e) 20 msec, (f) 23 msec, (g) 200 msec for the right column. 
The figure (d) and (h) shows the phase profile corresponding 
to Fig. (c) and (g), respectively. 



very close to each other in the central hole, they never 
overlap. 

As fl is further increased, all vortices generated from 
the condensate surface are absorbed into the central den- 
sity hole as shown in the right column of Fig. 7, where the 
central hole is composed of 24 singly-quantized vortices 
packed together. The packing is possible for high rota- 
tion frequencies (fl > uj±) because the centrifugal force 
decreases the central condensate density, and thus pack- 
ing the vortices together costs less energy. The minimum 
of Eq. (|14|l determines the radius of the ring condensate 
in Fig. 7(g) as i? = a/to(S12 — uj'j_)/k. Here we call a set 
of vortices such as that in Fig. 7(g) " a giant vortex" 
to indicate that a number of phase defects are contained 
in a single hole. Figures 7(h) also shows that the phase 
singularities do not completely overlap. 



13 



A voitex lattice 
£1=1.801 



A circular vortex 
aiTay with a hole 



A giant vortex 



n=2.5Ki 



n=3.0ai| 



® lolo 



n=2.o% 




n=2.75% £1=3.3% 



olo 



2.0 



3.0 



D/co 



FIG. 8: The equilibrium vortex states obtained in the numer- 
ical simulation for C = 250. 



The obtained vortex states for the coupling coefficient 
C — 250 are summarized in Fig. 8 as a function of a 
rotation frequency f2. For small f2, the equilibrium state 
is the vortex lattice state. As VI increases, the vortices 
begin to merge in the central region. We find that two 
vortices merge at the center at ~ 2.2w_l. This is the 
onset of the appearance of a new vortex state consisting 
of a density hole and a circular array of vortices. Further 
increase of il (above 3.0) stabilized the giant vortex. Re- 
cently, Kavoulakis and Baym studied the vortex states in 
a fast rotating regime under the Thomas-Fermi approxi- 
mation |25j, which is effective only for the large interac- 
tion energy C > 5000. It is difficult to make the numer- 
ical simulation for such a strong interaction because the 
high rotation will generates an extremely large number of 
vortices, whose numerical description would need a very 
high spatial resolution. It is interesting to investigate 
the phase diagram of rich vortex states in a fast rotating 
BEC; Kavoulakis and Baym proposed a phase diagram 
with triple points between three distinct vortex states in 
the interaction strength versus rotation frequency plane 

m 

Experimentally, this giant vortex has been created by 
Engels et al. [t^ . The authors removed atoms of a rotat- 
ing condensate from the central region, which resulted in 
the further increase of the angular momentum per parti- 
cle. Eventually, about 60 vortices merged into a density 
hole, which had the long lifetime of many seconds. The 
combined harmonic-plus-quartic potential was created in 
ENS group recently by superimposing a blue detuned 
laser with the Gaussian profile |77J. They observed that, 
as the rotation frequency approached the frequency of 



the harmonic potential, the clear image of a vortex lat- 
tice disappeared, which suggested a transition into a new 
vortex phase. 

V. CONCLUSIONS 

In this review article, we discussed the dynamics of 
quantized vortices in superfluid helium and atomic-gas 
BECs. Analogous topological defects appear in almost 
all fields of physics, including superconductors, magnetic 
materials, and cosmology. In each such field, the defects 
play important roles in the physics. However, unlike most 
systems in which the physics of topological defects usu- 
ally involve impurities, the vortices in superfluid helium 
and atomic BECs exist in very pure systems in which the 
proper physics of the topological defects themselves are 
relevant. This is the first reason why we are interested 
in vortices of these systems. The second reason is that 
quantum mechanical constraints make the vortices more 
stable and definite than those in classical media. 

It has been about 40 years since a quantized vortex 
was discovered in superfluid helium. As described in 
this article, research on superfluid turbulence is about 
to reach a new stage regarding the similarity between su- 
perfluid turbulence and classical turbulence, a stage in 
which the nature of quantized vortices is crucial. In ad- 
dition, the more recent realization of atomic-gas BECs 
have allowed us to visualize and control the dynamics of 
vortices, such that the quantum Hall state with a very 
high density of vortices may soon be reached [73, l7lLl72ll . 
Unconventional vortices in multicomponent BECs j23| is 
another area that will likely produce new and interesting 
results. Given the recent history of rapid developments in 
these fields, even more unexpected phenomena are likely 
to emerge in the future. 
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